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The Universe may harbor relics of the post-inflationary epoch in the form of a network of self- 
ordered scalar fields. Such fossils, while consistent with current cosmological data at trace levels, 
may leave too weak an imprint on the cosmic microwave background and the large-scale distribution 
of matter to allow for direct detection. The non-Gaussian statistics of the density perturbations 
induced by these fields, however, permit a direct means to probe for these relics. Here we calculate 
the bispectrum that arises in models of self-ordered scalar fields. We find a compact analytic 
expression for the bispectrum, evaluate it numerically, and provide a simple approximation that 
may be useful for data analysis. The bispectrum is largest for triangles that are aligned (have 
edges fci ~ 2&2 — 2&3) as opposed to the local-model bispectrum, which peaks for squeezed triangles 
(fei ~ ki ;$> k$), and the equilateral bispectrum, which peaks at k\ ~ &2 — k$. We estimate that this 
non-Gaussianity should be detectable by the Planck satellite if the contribution from self-ordering 
scalar fields to primordial perturbations is near the current upper limit. 

PACS numbers: 98.80.-k, 98.80.Cq, 11.30.Fs 



I. INTRODUCTION 

A wealth of precise cosmological data are in good 
agreement with the predictions of the simplest single-field 
slow-roll (SFSR) inflationary models [l]. Still, no theo- 
rist considers these as anything more than toy models. 
Realistic models must surely be more complicated, and 
they generically predict that there should arise, at some 
point, observable phenomena that depart from the pre- 
dictions of SFSR inflation. Some possible directions for 
physics beyond the SFSR approximation include multi- 
field models fei 3| and inflaton models with non-standard 
kinetic terms [4 1 . There has also been investigation of the 
consequences of topological defects [5| produced toward 
the end of or after inflation [6| . 

If inflation was followed by a transition associated 
with the breaking of a global O(N) symmetry, then self- 
ordering scalar fields (SOSFs) are another possibly ob- 
servable ear ly-Uni verse relic, even if there are no topo- 
logical defects (i.e., if N > 4). Here, the alignment 
of the scalar field as the Universe expands gives rise 
to a scale-invariant spectrum of isocurvature perturba- 
tions, without topological defects 0- Sample variance on 
the current data limit these perturbations to contribute 
no more than ~ 10% of large-angle cosmic-microwave- 
background (CMB) anisotropy power SOSF mod- 
els are parametrized simply by the number N of scalar 
fields and the vacuum expectation value v. The CMB 
constraint implies (v/N 1 / 4 ) < 5 x 10 15 GeV, as we ex- 
plain below. At this low amplitude, it is unlikely that 
any surviving relics leave a distinct imprint on the CMB 
power spectrum [ll| . 

In recent years, non-Gaussianity has been developed as 
a novel tool to investigate beyond-SFSR physics [l2.ll3j. 
SFSR models do not predict that primordial pertur- 



bations should be Gaussian, but the departures from 
Gaussianity that they predict are unobservably small 
[Ij-flQ]- Multi- field models 3, such as curvaton mod- 
els Q, string-inspired DBI 0, LL7| models, and models 
with features in the inflaton potential [HI, [H[ can all 
produce larger, and possibly observable, deviations from 
non-Gaussianity. For example, the detailed shape (trian- 
gle dependence) of the bispectrum may also help distin- 
guish these different scenarios. The "local-model" bis- 
pectrum, like that which arises in curvaton and multi- 
field models, has a very different shape dependence than 
"equilateral-model" bispectra, like those in DBI mod- 
els. Non-Gaussianity can be sought in the CMB [l9j . 
large-scale structure (LSS) [2(|, and the abundances and 
properties of gravitationally-bound objects plj or voids 
[22| . Biasing may significantly amplify the effects of non- 
Gaussianity [23| in the galaxy distribution. 

The energy-density perturbations in self-ordering 
scalar fields are quadratic in the scalar-field perturba- 
tion, which may itself be approximated as a Gaussian 
field. The density perturbations induced b y S OSFs are 
thus expected to be highly non-Gaussian 0,0, H^], even 
in the absence of topological defects. It is thus plausi- 
ble that the non-Gaussianity induced by SOSFs might 
be detectable, even if they provide only a secondary con- 
tribution to primordial perturbations. 

In this paper, we perform the first calculation of the 
full shape (triangle) dependence of the bispectrum from 
SOSFs. We follow the formalism for non-Gaussianity 
developed in Ref. [24j. We find considerably simplified 
formulas for the bispectrum, evaluate them numerically, 
and find a simple approximation to aid in data-analysis 
efforts. We estimate the current non-Gaussianity con- 
straint to the model parameter space and find it to be 
comparable to that from the upper limit to isocurvature 
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perturbations from CMB fluctuations. 

The plan of this paper is as follows: In Section [Til we 
define the model, write the scalar-field equations of mo- 
tion, show that the dynamics are those of a nonlinear- 
sigma model, and introduce the large- -/V scaling limit for 
the nonlinear sigma model. In Section IIII1 we write the 
relation between the matter-density perturbation and the 
scalar-field perturbation. In Section IIV| we derive the 
power spectrum for density and curvature perturbations, 
discuss the normalization, and derive current constraints 
to the v-N parameter space from upper limits to the 
SOSF contribution to CMB fluctuations. In Section fVl 
we discuss the calculation of the bispectrum, the cen- 
tral focus of this paper. We present a simplified ver- 
sion, our Eq. (|26[) . of the matter-bispectrum expression 
in Ref. [24j , evaluate it numerically, and provide a sim- 
ple analytic approximation for the results. We write the 
bispectra for matter and curvature perturbations, define 
a non-Gaussianity parameter for the model, and esti- 
mate the current constraint to from the CMB. Section 
IVII presents the matter bispectrum for modes that en- 
tered the horizon during radiation domination, those rel- 
evant for galaxy surveys. The central results of the paper 
are Eq. ([55)1 for the curvature bispectrum; Eq. (1521 which 
defines /°j in terms of the SOSF model parameters v and 
N; Eq. (|28[) which approximates the bispectrum func- 
tion 53(^1,^2,^3); and Eqs. (|34]) and (|35|) which present 
the matter bispectrum in a form useful for galaxy sur- 
veys. We make concluding remarks in Section IVIII An 
Appendix contains some calculational details and useful 
approximations. 



II. SCALAR-FIELD DYNAMICS 

Self-ordering scalar fields are described by an TV- 
component scalar field with an O(N) symmetry that is 
spontaneously broken to 0(N — l). 1 After symmetry 
breaking, the scalar field lies in different places in its 
S 1 ^^ 1 vacuum manifold in different causally disconnected 
regions of the Universe. As the Universe expands and 
these previously causally-disconnected regions come into 
causal contact, field gradients tend to align the scalar 
field. The rate of alignment for these fields is limited only 
by causality, and so the fields become aligned within a 
few Hubble times after horizon crossing. Still, as the 
Universe expands, there are continually new causally- 
disconnected regions, on ever larger scales, that enter 
the horizon. The result is thus a continual scale-invariant 
generation of new scalar-field perturbations. In this Sec- 
tion, we describe the scalar-field dynamics; the following 
Section then describes how the gradient energy density 
in these scalar fields induce perturbations to the matter 



1 We assume that the issues about global symmetries raised in 
Refs. 26] are somehow solved 12711. 



density. 

The starting point is an iV-component scalar field $ = 
((f) 1 , (f> 2 , • • • , <fi ), with (f> a real, with Lagrangian density, 



(1) 



where A is the dimensionless self-coupling of and v 
is the magnitude of the vacuum expectation value (vev) 
in the true vacuum. At temperatures T <C X 1 ^ 4 v, the 
O(N) symmetry of the Lagrangian is spontaneously bro- 
ken, and the field is thereafter restricted to the S 1 ^^ 1 
vacuum manifold. The dynamics is thus effectively that 
of N — 1 massless Nambu-Goldstone modes which we 
describe in terms of the N fields <fi a with the effective 
Lagrangian density, 



C = -(V M $) • + A(|$| 2 - v 2 ), 



(2) 



where A is a Lagrange multiplier that enforces the con- 
v 2 . The resulting equations of motion are 



straint |$| 



a "(x,ry) + W(x,tj) 



V 2 + (V M $) • 



0°(x,»j)=O, 



(3) 



where the primes denote derivatives with respect to con- 
formal time 77, and % = a' /a in terms of the Friedmann- 
Robertson- Walker scale factor a(rj). Also, V 2 is here a 
spatial Laplacian in comoving coordinates. Eq. ([3]) repre- 
sents the non- linear sigma model (NLSM from now on), 
that describes the evolution of the scalar field after spon- 
taneous symmetry breaking. 

In the large- TV limit, the field components become in- 
dependent of each other (up to corrections of order A -1 ). 
We thus replace the bilinear term in the equation of mo- 
tion by an ensemble average, 

(V M $) • (V"l) = N ((V^°)(V^°)) = T( v ) , (4) 

where there is no sum on a in the second equality, and 
in the last equality we have made the usual ergodic as- 
sumption, replacing the ensemble average by a spatial 
average. 

The only timescale in the problem is that set by the 
(comoving) horizon 'Hr x oc 77, so by dimensional con- 
siderations T oc H 2 , and T(r)) = T /r] 2 , with T Q > 0. 
We then replace the non-linear term in the NLSM equa- 
tion of motion, Eq. ([3]), by this expectation value and in 
this way linearize the equations of motion. Introducing 
a = dloga/dlogi] and Fourier transforming the spatial 
dependence of the equations, 



we obtain 



6 a CM) 



2a 



c? 3 x0 a (x, 77) e 



-Hk-x 



+ \k z - 



T 
v 2 rj 2 



^ = 0, 



(5) 



(6) 
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with a = 1 for a radiation-dominated Universe and a = 2 
for a matter-dominated Universe. For constant a, the 
solution to Eq. ([6]) that is finite as 77 — > 0, is </> a (k, rf] = 
w a (k)/(fc?7), with /(a;) = x 1 / 2 "" J„(a;), and J„(x) is a 
Bessel function. Here, e°(k) is the amplitude of mode k, 
and v is fixed by v 2 = (1/2 - a) 2 + (T /w 2 ). 

In the large- N limit, the statistical distribution of each 
field component approaches a Gaussian distribution with 
mean (</> a (x, rj)) = and variance (</> a (x, rj)4> h {x., rjf) = 
(v 2 /N)8 a t- The initial field component ^ a (x, 77 = 0) takes 
on a random value at each point in space. We thus take 
the {e°(k)} to be Gaussian random variables with mean 
(e a (k)) = and variance, 

(e a (k)e\k')) = (2vr) 3 |kr|^Mk + k'), (7) 

where <$o(k) is the Dirac delta function, and A is a nor- 
malization constant to be determined below. The power- 
law dependence on k is taken since the initial conditions 
are scale-free. 

The power-law index n in Eq. (JTJ) is fixed by the con- 
dition that (0 a (x,77)/(x,7?)) = (v 2 /N)5 ab for all 77: 



6 Q (x,7 7 )/(x, ?7 )) = S ab v 



d 3 k f d 3 k' 



(2tt) 3 J (2tt) 3 
x (e Q (k)e Q (kO)/(fc77)/(fc'77) 

^TTV 2 S ab 



{2n) 3 AN 



dfcr-"/ 2 (fcr/).(8) 



We see that n — 3 gives a result that is independent of 
time, and so we choose n — 3 hereafter. 

Just after symmetry breaking, at conformal time 77*, 
the field correlation is then 

(<T(k,r?»)/(k',77»)) ex f 2 (x*)(e a (k)e b (k>)) 

oc ^M-^-Hnik + k'). 

(9) 

Since the initial values a (x, 77*) are uncorrelated on 
scales k <C 777 1 , we set ^ = a + 1, so that the initial 
field is described by a white-noise power spectrum. This 
then fixes (T /tj 2 ) = 3a + (3/4). 

We now return to Eq. (JSJ to fix the normalization con- 
stant A. From 



<0 Q (X,77)/(X,77)) = 



ab 



2tt 2 AN 



dxx- 2a J 2 a+1 {x) 



V 2 Sab 

N 



we find 



.4 



T(a) 



87r 2 r(2a + 3/2)r(a+l/2) 



(10) 



(11) 



For a — 2 (matter domination), A = 16/28357T 3 = 
1.82 x 10 -4 , and for a = 1 (radiation domination), 
A = 2/15tt 3 =4.3 x 1CT 3 . 



III. MATTER-DENSITY PERTURBATIONS 

Although the scalar field will initially take on differ- 
ent values in different causally-disconnected regions, the 
curvature perturbation is initially zero. The scalar-field 
gradient-energy perturbation that arises as previously 
causally-disconnected regions come into causal contact is 
then compensated by a perturbation in the matter den- 
sity [13,[2i|. 

In this Section, we calculate the time evolution of the 
matter perturbation. The action of the scalar field oc- 
curs primarily within a few Hubble times after a partic- 
ular Fourier mode k enters the horizon. The subsequent 
evolution of the mode is then governed by gravitational 
infall as if it were a primordial perturbation; i.e., the 
perturbation amplitude grows only logarithmically dur- 
ing radiation domination, and then grows with the scale 
factor during matter domination. Our strategy here will 
be to evaluate the matter-perturbation amplitude several 
Hubble times after horizon crossing, a calculation that is 
relatively straightforward. Strictly speaking, our calcula- 
tion applies only to modes that enter the horizon during 
matter domination, but we argue below that our ulti- 
mate results for the bispectrum should also be roughly 
valid for the smaller-scale modes that enter the horizon 
during radiation domination, those relevant for galaxy 
surveys. 

As described in Ref. 24(, the scalar- field align- 



ment involves density perturbations that then lead to 
gravitational-potential perturbations which in turn in- 
duce the perturbations to the matter density that are our 
ultimate interest. Following Ref. (24 |. the matter-density 
perturbation induced by the scalar field for modes that 
enter the horizon during matter domination is 



2ttG 



?f / drfdiToiM): 



<5(x,77) 



where G is Newton's constant, and 

t 0i = (dormr) ■ 



(12) 



(13) 

is the Oi-component of the stress-energy tensor of the 
multicomponent scalar field. The integral in Eq. (fT2|) 
approaches a constant for 77 ^> iew/k — i.e., within a few 
Hubble times after horizon crossing. The subsequent rj 2 
evolution in the prefactor is then simply the S oc a oc tj 2 
linear-theory growth of the perturbation amplitude in a 
matter-dominated Universe. 

Using G = 1/M|i, where M P1 = 1.22 x 10 19 GeV is the 
Planck mass, and defining 



C 



2n f v 



5 \Mpi 



(14) 



the Fourier transform of the density perturbation is 
£(k,77) = -C77 2 |^ e a (q) e a (k-q)|k-q|(k.q) 



(2tt)3 

dr/'(|k-q|r)/(gr), 



(15) 
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where f'(y) = df/dy. The crucial qualitative feature is 
that 5(k, 77) is quadratic in powers of e a (k). And since 
e a (k) is a nearly Gaussian field, the density field 6(x) will 
be highly non-Gaussian. 



IV. THE POWER SPECTRUM 

The power spectrum P a {k) for matter-density pertur- 
bations induced by the scalar field is defined by 

(S(k)S(k')) = (2tt) 3 S D (k + k') P° (k), (16) 

where the angle brackets denote an average over all re- 
alizations of the random field 6(k). The calculation of 
the power spectrum is lengthy but straightforward; de- 
tails are provided in the Appendix. The result, given in 
Eq. (j40|) . can be re- written, 



P°{k,rj) 



C*V k 
(2n) 2 A 2 N J 



dvv i 



dll(v,b)l [l(v,b)vl+I(b,v)(l-vl)] , 



where b = vl 



2v L and 



J(a, b) 



ds 



f(as)f(bs) 

a 3/2 6 l/2 



(17) 



(18) 



Strictly speaking, the upper limit in this integral is krj. 
However, here we will restrict our attention to modes 
that have evolved well within the horizon, i;?; > 1, and 
so we take the upper limit of the integral in Eq. (|T5)) to be 
infinity. In this case, the integral I(a, b) is antisymmetric 
in its arguments, and the power spectrum can be written 



P a (k, V ) 



CV k 
A 2 N 



.92. 



(19) 



where 



.92 



d 3 v 



[I( W ,|z-v|)] 2 (z-v) [2(z.v)-l] , (20) 



and z is a unit vector. Details on the evaluation of X 
are given in the Appendix. For a — 2 (matter domi- 
nation), the integral evaluates to .g 2 = 3.3 x 10- 7 and 
for a — 1 (radiation domination) it is 52 = 2.1 x 10~ 4 . 
Note that the ratios $2 /A 2 that appear in Eq. (fl~9|) are 
approximately 10 and 11, respectively, for a = 2, 1, im- 
plying that the amplitude of the matter perturbation in- 
duced by the unwinding of the scalar field is the same, to 
0(10%), for modes that enter the horizon during matter 
and radiation domination. 



CMB temperature power spectrum at multipole moment 
I = 10. 

On subhorizon scales during matter domination, the 
curvature perturbation C( x ) is related to the gravita- 
tional potential $(x) by £(x) = (5/3)$(x). The grav- 
itational potential is related to the density perturbation 
through the Poisson equation, V 2 $ = AirGa 2 p8, where 
p is the mean density. In Fourier space, the curvature 
perturbation £(k) is thus related to the matter-density 
perturbation 5(k) by 



C(k) 



5 (aB 
2 \~k 



*(k) 



(21) 



where we have used the Friedmann equation H 2 = 
87rGp/3, and H = (da/dt)/a is the expansion rate. The 
amplitude of the curvature power spectrum due to the 
SOSF is therefore, 



A 2 



k 3 k 3 
2^ = 2^ 

^ 4 

92 



pi 



A 2 N 



•5 
2 

: 80 



aH 
~k~ 



P(k) 



M, 



pi 



1 

N 



(22) 



where we have used (aHaf) = 2 during matter domina- 
tion. The next step is then to determine the relation 
between the curvature-power-spectrum amplitude A 2 ^ 
and the temperature-fluctuation amplitude. This is a 
notoriously difficult calculation, but to get an estimate, 
we use Fig. 4 in Ref. [2{|, which shows that the large- 
angle temperature fluctuation AT in a SOSF model is 
G sw ~ 10 times greater than it would be in an adiabatic 
model with the same matter-power-spectrum normaliza- 
tion on large scales. 2 Current CMB measurements in- 
dicate a curvature power spectrum A-r. ~ 5 x 10~ 5 , if 
primordial perturbations are adiabatic. If the SOSF pro- 
vides G sw times more AT for fixed curvature, and if they 
make a fractional contribution p,y to the large-angle tem- 
perature variance, then A 2 ^ = (p (T /G 2 w )A 2 ^. We thus 
obtain 



AT 1 / 4 



P*A 2 A 2 , 
8G 2 W 52 



1/4 



Mi 



pi 



2000 ' 



(23) 



where the numerical result is obtained by taking p a = 0.1 
and G sw = 10. The numerical upper limit in Eq. (1231) is 
in good agreement with limits obtained from simulations 



A. Normalization of the power spectrum 

We now estimate the constraints to the v-N parameter 
space from the empirical constraint that the SOSF pro- 
vide no more than a fraction p a ~ 0.1 to G; = io @, the 



2 The factor of 10 is a bit larger than the factor of 6 one might 
attribute due to the difference (1/3 versus 2) for the Sachs-Wolfe 
amplitude for adiabatic and isocurvature perturbations. The ad- 
ditional AT may be due in part to the vector and tensor pertur- 
bations that are also excited in SOSF models. 
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V. THE BISPECTRUM 
A. The calculation 

The calculation of the bispectrum proceeds analo- 
gously. The matter bispectrum B(k\,k 2 ,k 3 ) is defined 
by 

(3(^)5(^)6^)) = (27T) 3 S D (-k 1 + k 2 + k 3 )B(fc 1 , fc 2) k 3 ). 

Although the definition of the bispectrum is nominally in 
terms of the vector quantities k^, the triangle constraint 
ki + k 2 + k 3 =0 imposed by the Dirac delta function, 
as well as statistical isotropy, imply that the bispectrum 
is most generally a function of the magnitudes ki of the 
three sides of the triangle. Again, some details of the 
calculation are provided in the Appendix. The result is 



C 3 rf 

B(k 1 ,k 2 ,k 3 ) = 3 2 ff 3 (fci,fc2,fc3) 



(25) 



where 

g3{ki,k 2 ,k 3 ) = 



d 3 v 



H(u + v,v) 



{2n) 3 ' 

xH(y,z- v)iJ(z- v,u + v), (26) 



with 



H(a,b) = l{a,b){b 2 -a 2 ). 



(27) 



We have chosen k\ in Eq. (|2"6")l to be in the z direction, 
without loss of generality, and we have then defined u = 
k 2 /ki. Note that H(a,b) — H(a,b) — H(b 7 a); i.e., it is a 
function only of the magnitudes of its arguments, and it is 
symmetric in its arguments. Note further that H (a, b) < 
0, and thus g 3 (ki, k 2 , k 3 ) < 0. The function g 3 (ki, k 2 , k 3 ) 
depends only on the shape of the triangle, not on its 
overall size — i.e., g 3 (ki,k 2 ,k 3 ) = 53(1, fo/^i, ^3/^1) — a 
consequence of the scale- invariance of SOSFs. We have 
checked that Eq. (|26]l is equivalent to, although far sim- 
pler, than Eq. (59) in Ref. ;24j. Given the symmetry 
of H(a, b) in its arguments, it is simple to check that 
g 3 (w — z — u) = <?(u), as it should (given that the three 
sides of the triangle should add asz + u + w = 0). If 
we set the third side to have length w = k 3 /ki, then 
cos 6 = u • z = (w 2 — 1 — u 2 )/(2u). If we choose 
k\ > k 2 > k 3 , then cos < — (2u) _1 . 

We have calculated g 3 {k\,k 2 ,k 3 ) numerically, and the 
result is shown in Fig. [TJ We note (prefacing the dis- 
cussion below) that the quantity, — 33(1, x 2 , x 3 ), with 
x 2 = k 2 jk\ and x 3 = k 3 /k\, that we plot is the 
same (up to some normalization factor) as the quantity 
F(q, x 2 ,x 3 )x 2 x 3 plotted in Figs. 1 and 2 in Ref. [32| which 
show, respectively, the bispectra for the local-model and 
equilateral model. Those figures show that the local- 
model bispectrum peaks sharply for "squeezed" triangles 
(fci ~ k 2 3> k 3 ) and that the equilateral-model bispec- 
trum peaks at equilateral triangles (fci ~ k 2 ~ fc 3 ). Our 




-83/A 3 



FIG. 1: The function —33(^1, fe, k^), taking ki — 1, for modes 
that enter the horizon during matter domination. The figure 
looks virtually identical for modes that enter the horizon dur- 
ing radiation domination. 



Fig.[T]shows that the SOSF bispectrum is, however, quite 
different. It is nonzero for equilateral triangles, goes to 
zero in the squeezed limit, and it peaks for "aligned" 
triangles, k\ ~ 2k 2 ~ 2A: 3 . 

To aid in data-analysis efforts, we have found that the 
following approximation reproduces the numerical results 
for <? 3 (fci, fc 2 , k 3 ) to within a few percent: 



g 3 (ki,k 2 ,k 3 ) = - 



A 3 / 






262 - 127^ 


143 \ 




947^ 


-"*>(£) 


k\ 





893 



(28) 



where we take k\>k 2 > k 3 in this expression. 



B. Curvature Bispectrum 

To compare with results for other models, and for com- 
parison with CMB constraints, we next calulate the cur- 
vature bispectrum F{k\, k 2 , k 3 ), defined by 

(C(ki)C(k 2 )C(k 3 )) = (27r) 3 fe(ki+k 2 + k 3 )F(fci,fc 2 ,fc 3 ). 



Using the relations above, we can write C 2 
(-K 2 /50)(A 2 NA 2 Za /g 2 ) and then find, 



(29) 



F ff (fci,fc 2 ,fc 3 ) = 



2V2ir 3 A 3 na g 3 {k x .. 



k 2 ,k 3 ) 



g 3 2 /2 N^ 



h 2 h 2 h 2 
/v 1 fi, 2 ft, 3 



(30) 



A few observations: (1) It is only the amplitude, not the 
shape, that depends on the symmetry-breaking scale v. 
(2) The bispectrum decreases, for fixed An a , as N^ 1 / 2 
with increasing TV, again reflecting that the model should 
become increasingly Gaussian with more fields, a conse- 
quence of the central-limit theorem. (3) The scaling with 
Ajz is oc A^j, as opposed to the A^ scaling of the local- 
model bispectrum. In words, the non-Gaussianity is of 
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order unity, a consequence of the fact that the density 
perturbation is the square of a Gaussian random field 
[cf., Eq. (|12p ]. rather than something very small, as in 
inflationary models. 

We now put the curvature bispectrum in a slightly 
more familiar form by defining a non-Gaussianity pa- 
rameter /°[ for the model. To do so, we recall that the 
local-model prediction for the curvature bispectrum is, 



F local (k 1: k 2 ,k 3 ) 



2^( 27 r 2 ) 2 A^ 
1 1 



local 
1 



1.31.3 
1*2 



+ 



"-2 "-3 



+ 



1.3 1,3 
K 1 K 3 



,(31) 



where / n ocal is the local-model non-Gaussianity parame- 
ter, defined by writing the curvature as £(x) = Cg( x ) + 
(3/5)/^ cal [(Cs) 2 (x) - <(Cs) 2 (x))] in terms of a Gaussian 
field C ff (x). 

We now define the non-Gaussianity parameter /°j for 
SOSFs by equating the local-model and SOSF bispectra 
for equilateral triangles; i.e., equating Eqs. (pTLj) and (1501) . 
we define, 



/nl 



5^ /2 53 (l,l,l) 



40 G 



-3 



18\/2^iVi/2 32 3 /2 ATCG 3 w 

, 3 / 2 / N\ ~ 1 I 2 
6 

„ / n \ 

3G 



-f-3 

' sw 1 5 x 10 15 GeV 



(32) 



With this f°i, the curvature bispectrum can then be writ- 
ten 



1 R 



2,-1 ,4 /nl g3(kl,k 2 ,k 3 ) 

k\k\k\ 5 3(1,1,1) 



(33) 



Similarly, the gravitational-potential bispectrum can be 
obtained by multiplying this expression by 5 /3 and then 
replacing A|, by A| = (3/5) 2 A^. Note that is mani- 
festly positive, unlike f n \ for the local or equilateral mod- 
els, which may take on either sign. 



C. Estimate of CMB Constraints 

As indicated in the Introduction, the bispectrum can 
be probed with the CMB, large-scale structure, and the 
abundances of objects. The strongest current constraints 
to the local-model bispectrum come from the CMB [33j . 
followed closely by galaxy-clustering constraints. Given 
that SOSFs produce a larger temperature-fluctuation 
amplitude for a given density-perturbation amplitude, we 
surmise that the CMB will provide stronger constraints 
to SOSF non-Gaussianity than galaxy clustering. We 
thus now estimate a constraint to from the CMB. 

Before doing so, we first caution that Eq. (|33p is de- 
rived for the curvature perturbation only for modes once 
they are well within the horizon. It is thus not, strictly 



speaking, appropriate for CMB modes I < 100. Still, the 
shape dependence of the bispectrum, and its amplitude 
relative to the curvature-perturbation amplitude, arises 
primarily from the quadratic dependence of the density 
perturbation on the scalar-field perturbation as encoded 
in Eq. (|15p . The shape dependence of the bispectrum we 
calculate should thus be at least roughly correct even for 
I < 100. 

CMB constraints to / n i are typically applied assum- 
ing that the density perturbations are adiabatic, which 
implies a certain relation, (AT/T) ~ — £/5, for the large- 
angle temperature fluctuation. In our case, though, there 
is roughly G sw ~ 10 times more AT for a given £ than 
in adiabatic models. If so, and if all of these tempera- 
ture fluctuations are due to scalar perturbations (rather 
than vector and/or tensor modes), then the implied CMB 
bispectrum should be roughly G^ w times larger. Simu- 
lations show, though, that only a fraction / s ~ 0.5 of 
the SOSF temperature-fluctuation power is due to scalar 
modes, the rest coming from vector and tensor perturba- 
tions [1^, |3(| • The implied CMB bispectrum should thus 
scale with / s 3 ^ 2 . Combining these scalings with /°j oc G^ 
[see Eq. ([3"2|] , the G sw dependence of the CMB bispec- 
trum drops out. We can therefore apply CMB constraints 
to /^j by identifying the f n \ constraints obtained from 

the CMB for adiabatic perturbations with Gg W / s 3/ ' 2 / n T 1 . 
And one final caveat: We disregard the differences in the 
temperature power spectra from SOSFs and adiabatic 
perturbations. 

Keeping these multiple caveats in mind, we proceed 
with our very rough estimate by noting that the WMAP- 
7 95% C.L. limit to /* qul \ the non-Gaussianity parameter 

for the equilateral model, is —211 < /"j 51 " 1 < 266 [3J|; this 
bispectrum is maximized for equilateral triangles. On 
the other hand, the SOSF bispectrum is maximized for 
aligned triangles and is zero for squeezed triangles. We 
thus conclude that the constraint to Gg W / s 3//2 will be 
stronger than that to /^ qul1 , but it is not clear — given the 
different weightings to squeezed and aligned triangles — 
how it will compare with that to / n ocal . Applying these 



rough arguments to Eq. (|32l) . with f s ~ 0.5, then we 
estimate a non-Gaussianity parameter in excess of the 
predicted threshold f n \ ~ 7 for detection by Planck [35j . 
For now, however, we simply estimate conservatively that 

\o 3 sw fs /2 a < 200. 



VI. THE BISPECTRUM FOR GALAXY 
SURVEYS 

We have carried out our calculations in the regime 
where analytic progress is most easily made — i.e., modes 
that have entered the horizon during matter domination 
and only after those modes have evolved well within the 
horizon. Strictly speaking, therefore, our calculations ap- 
ply only to galaxy surveys on very large scales — those 
generally larger than extant surveys cover — and possibly 



7 



to mid-scale regimes in the CMB. 

Still, our results for the bispectrum can be easily 
adapted to obtain roughly the bispectrum for smaller- 
scale modes, those that enter the horizon during radia- 
tion domination and those relevant for galaxy surveys. 
The calculation for the evolution of these modes is al- 
tered by three effects: (1) The index v = 1 + a for the 
Bessel functions in the scalar-field dynamics ( Section HI]) 
is different. However, we have shown that this has no 
more than an 0(10%) effect on the power-spectrum and 
bispectrum. This is simply because the correlations be- 
tween different Fourier modes of the scalar-field energy 
density are imprinted, through Eq. (fT5|) , by the depen- 
dence of those Fourier modes on the scalar-field perturba- 
tions. This is the same for modes that enter the horizon 
during matter and radiation domination. The other two 
effects are (2) a slightly different amplitude for the matter 
perturbation, relative to the scalar-field energy-density 
perturbation, for modes that enter the horizon during 
radiation domination [Eq. (32) for RD in Ref. 24], as op- 
posed to Eq. (29) in the same reference, our Eq. (|12l) . for 
MD]; and (3) the usual linear-theory growth of primordial 
isocurvature perturbations through radiation domination 
and the transition to matter domination. These latter 
two effects amount to a calculation of the transfer func- 
tion T(k) for the matter power spectrum in SOSF mod- 
els, which can be accomplished either with simulations 
[cf., Ref. [29|, [3l[] or approximately with standard linear- 



theory calculations with primordial isocurvature fluctu- 
ations. Again, however, although the calculation of the 
evolution of the amplitudes of the small-scale density- 
perturbation Fourier modes will be far more complicated 
than the larger-scale modes we have focused upon, the 
correlations between those modes will be, at the 0(10%) 
level, the same as those we have calculated for larger- 
scale modes. 

More precisely, all we need to do is replace the density 
fields S(k) in Sections JV] and by 5(k)T(k), where 
T(fc) is the SOSF transfer function. The matter power 
spectrum P a {k) due to SOSFs is then obtained from that 
in Eq. (IT9l by multiplying by |T(fc)| 2 , and the matter 
bispectrum is obtained by multiplying that in Eq. (|25|) by 
T {k\)T (k 2 )T (k 3 ) . We can then write the normalization 
constant (Cij 2 /A) 3 in Eq. (|2"5|) in terms of the (processed) 
matter power spectra using Eq. (TT9"]) to obtain the matter 
bispectrum, 



B(k u k 2 ,k 3 ) 



g 3 (fci,fc 2 ,fc 3 ) 



P CT (fci)P CT (fc 2 )P CT (fc 3 ) 
k\k 2 k 3 



1/2 



valid for galaxy-survey scales. Here P a {k) is the pro- 
cessed power spectrum due to SOSFs; i.e., it includes the 
transfer function. Using Eq. (1321) . this can be re- written 
in terms of as, 



B(ki,k 2 ,h) = 



18V2tt fZAnG 3 sv , g 3 {k x ,k 2 , k 3 ) 



25/5 



■yp, 

G,. 



3/2 



53(1,1,1) 



P°(k 1 )P°(k 2 )P°(k 3 ) 



kik 2 k 3 



1/2 



SW 

1CT 



- 3 / 2 gafol,**,**) 
53(1,1,1) 



P g (fcl)P J (fc 2 )-P CT (fc 3 ) 

k\k 2 k 3 



-1 1/2 



(35) 



We leave further evaluation of this bispectrum, as well as 
assessment of current constraints, for future work. 



VII. DISCUSSION 

If some post-inflationary physics involves the sponta- 
neous breaking of an exact 0{N) symmetry with N > 4, 
then the ordering of these scalar fields may provide a sec- 
ondary contribution to primordial perturbations. Cur- 
rent constraints allow up to ~ 10% of the power in pri- 
mordial perturbations to be due to SOSFs. SOSF models 
are appealing from the theoretical perspective because 
they are simple, well-defined, and parametrized only by 
the symmetry-breaking scale v and number N of fields. 

In this paper we have calculated the matter and cur- 
vature bispectra induced by the ordering of such scalar 
fields. Given that the density perturbation is quadratic 



in the scalar-field perturbation, SOSF density perturba- 
tions are expected to be highly non-Gaussian, and if so, 
measurements of non-Gaussianity may provide the means 
to test these models. 

Here we have calculated analytically the bispectrum 
due to SOSFs and presented results in a way that 
should be easily accessible to those doing measurements 
with the CMB and large-scale structure. We find that 
the triangle-shape dependence of the bispectrum peaks 
for aligned triangles, unlike the local-model bispectrum, 
which is largest for squeezed triangles, and the equilat- 
eral bispectrum, which is largest for equilateral trian- 
gles. We have estimated a current upper limit to the 
non-Gaussianity parameter for the model and find 
that the implied constraints to the v-N SOSF parameter 
space are competitive with those from the upper limit to 
CMB temperature fluctuations. 

Finally, we have already argued above, in Section IV CI 
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that the correlation of modes will be similar for the large- 
scale modes as they enter the horizon, those relevant for 
large-angle CMB fluctuations. We therefore believe that 
rough constraints to the model can be derived from CMB 
measurements by assuming that the curvature bispec- 
trum we calculate is the primordial one. 

Clearly, there is room for further numerical work to 
test our assumptions and to make our predictions more 
precise. In the meantime, though, we believe that our an- 
alytic approximation captures the essential physics and 
that our bispectrum can be used in the meantime as a 
"working-horse" model to derive constraints, from non- 
Gaussianity measurements, to this interesting class of 
models for secondary contributions to primordial pertur- 
bations. 

Finally, we note that the model makes a number of 
other predictions. Given that density perturbations are 
actively generated as new modes come within the hori- 
zon, vector and tensor modes will be excited, and these 
may give rise to interesting polarization signals [36[ in 
the CMB and perhaps excite B modes 37J in the CMB 
that might be distinguished from those due to inflation 
[38j . There will also be a scale- invariant spectrum of pri- 
mordial gravitational waves produced [39| that can be 



sought in gravitational- wave observatories. 
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Appendix: Calculational Details 
A. The Power Spectrum 



From Eq. (|T5|) . the two-point correlator in Eq. (fl6|) is 

expressed in terms of the correlator of four e a (k)'s. Since 
e a (k)'s are gaussian distributed, we find via the Wick 
theorem, 



(5(k)5(k>)) = CV/^(^(qK(k-q) £ 6 (q') e fc (k'-q')) 

x|k-q||k'-q'|(k-q) (k' • q') J dr /'(|k - q|r)/(«r) / de/'(|k' - q!\OfWO 
= J ^0- [< e a (qK(k - q)> < £ b (qV(k' - q')) 

+ (e a (q)e b (q')) (^(^ - ^"(V - q')> + (^(^O^' - q')) (^(^ - q>V)>] 
x|k-q||k'-q'|(k-q) (k' ■ q') / dr /(«r)/'(|k - q|r) / dtftfZ)f(\V ~ q'|£) 



Using Eq. 0, we find 



and 



7VAV|k-q| 3 



(k • q)(k' • q') (e a (q)e h (k' - q')> (e a (k q) e b (q')) 

(27T) 6 



NA 2 q 3 \k-q\- 



(36) 



(k • q)(k' • q') (£ a (q)£ a (k - q)) {e b (<lY(k' - q')) 

= (k • q)(k' • q') J^ 3ql3 S ^b b S D (k)S D (k') = 0, (37) 



(k • q)(k' • q') (e a (q)e 6 ( q ')) (e a (k - q)e b (k' - q')> 
(2nf 



(k.q) 2 fe(q + q')fe(k + k'), (38) 



(k ■ q) (k ■ (k - q)) S D (q! - q + k)6 D (k + k') . (39) 
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Substituting these expressions into Eq. (|36|) . we obtain 



P°{k) = 



CV f d 3 q 



A 2 N J (2tt) 3 



dr 



f( q r)f'(\k- q \r) \ 2 
dT g3/ 2 | k _ q |i/2 ) !k ' ( f 

/(gr)/'(|k-q|r) 



g 3 / 2 |k — ql 1 / 2 



cV k r d 3 



X(v, |k - v|) (k • v) 2 X(v, |k - v|) + (k ■ v)(l - k • v)X(|k - v|, v) 



A 2 N J (2tt) 3 " 

where X(a, 6) is defined in Eq. (|18p . From here we can then introduce g 2 in Section ITVl 



(40) 



B. The Bispectrum 



The bispectrum is obtained by starting with, 



(6(^)6(^)6^)) = -C 3 V 6 J d3<ll ^ ) 2 f q3 (e a (qi) £ a (kl - qi) £ b (q 2 )e"(k 2 - q 2 )e c (q 3 K(k 3 - q 3 )) 
x |ki - qi||k 2 - q 2 ||k 3 - q 3 |(ki • qi) (k 2 • q 2 )(k 3 • q 3 ) 

x /rfn/(«m)/'(|ki -qi|n) / dr 2 /(g 2 r 2 )/'(|k 2 - q 2 |r 2 ) J dr 3 /(g 3 r 3 )/'(|k 3 - q 3 |r 3 ) . 

(41) 

The expectation value of the product of six e a (k)'s can be expanded with Wick contractions and, after some algebra, 
and using Eq. (O, we find 



(5(k 1 )5(k 2 )5(k 3 )) = -(2^) 3 <5 D (k 1 +k 2 + k. 



C 3 rf f d 3 v 



(k-v) l(v,b) 



A 3 N 2 J (2tt) 3 

x { [(u • v) l(v, 6 2 ) - (u • b 2 ) l{b 2 ,v) } [ ((k + u) • b) 1(6, fe) + ((k + u) • b 2 ) X(6 2 , 6) ] 

+ [(u • b X(6, 6i 2 ) - (u ■ b u ) 1(612, &) ] [ ((k + u) • v) X(«, 6 12 ) + ((k + u) • b 12 ) X(b 12 , v) } } , 

(42) 

where k = |ki |, k = ki/fc, u = k 2 //c, and we have defined b = k — v, b 2 = u + v, bi 2 = k + u v. Defining <? 3 as in 
Eq. ([23)) and performing a change of variables v — > (k — v) in the second term of the above integral, we then find 



ga(k 1 ,k2,k 3 ) 



^3 [(k • v)X(«, 6) + (k • b)X(6, «)] [(u • v)X(v, 6 2 ) - (u ■ b 2 )X(6 2l «)] 



[b • (k + u)X(6, 62) + b 2 • (k + u)X(6 2 , 6)] 



(43) 



If the upper limit of the integral defining I(a, b) is much greater than one (as expected for sub- horizon modes), such 
that X(a, b) = 1(6, a), we then obtain 



93 



d 3 v 



■ l{v, 6) l(v, 62) X(6, 62) (k • (v - b)) (u • (v + b 2 )) ((k + u) • (b - b 2 )) 



(44) 



Finally, simple algebraic re- arrangements in Eq. (144)) then results in the far simpler expression, in Eq. (|26p for 173 

d 3 v 



.93 



(2tt) : 



■ H(u + v, v) H(y, k-v) H(k - v, u + v). 



(45) 



with H(a,b) = (6 2 — a 2 )l(a, 6) < 0. This expression is indeed equivalent to Eq. (59) in Ref. [24|, althou gh i t is 
written in a much more compact and simpler way. Note that it is smaller by a factor (27r) 3 than that of Ref. [2J|, as 
a consequence of different conventions. 



C. Some Integrals and approximations 

Once modes are well inside the horizon, the upper limit 
for the integral in Eq. (fT5)l is large, and the integral can 



then be approximated by 
X(a, 6) e 



ds 



a 3 / 2 6 1 /2 



(46) 
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where f(x) = x 1 / 2_Q Ji+ Q , and a = 2 for matter domi- 
nation. The integral can be performed analytically; the 
result is 



^F(a 2 /b 2 ), for a < b, 



X(a, b) = 



where 



F(x) 



4" 



^F(b 2 /a 2 ), for a > b, 

X2-F1 (§, § - n;n + 2;x) 
r(n + 2)r(ra- §) 



2F1 (f, § - + l;a:) 
r(n + l)r(n-|) 



(47) 



(48) 



2-F1 (w, x; y; z) is the hypergeometric function, and n — 
1 + a is the order of the Bessel function. While straight- 
forward to evaluate numerically, this exact solution may 
be computationally expensive to evaluate repeatedly. We 
therefore use for our numerical work the approximation, 



-1 (b/a)"-l - f 1 

oBP' {b/a)»+l> 11 a °> 
I(a,6)~^ (49) 

1 (a/b)"-l jf h . 

. 96^ (a/b)« + l ' 11 ^ 



which provides good agreement with the exact results 
with k = 2.5. 
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